Scripts Matheus Sousa

					#105.1
#lendo os dados
esaligna<-read.table("esaligna.txt", head=TRUE, as.is=TRUE, sep="," )
#abrindo o arquivo
jpeg(filename = "Dispersao_105.1.jpg", width = 2000, height = 800, units = "px", pointsize = 12, quality = 100, bg = "white",  res = NA)
#plotando
par(tck=(0.03))
plot(dap~ht, data=esaligna, xlab="Altura Total (m)", ylab="Diametro a alatura do peito (cm)", main="Dispersao de DAP por altura total em E. Saligna ", cex=1.2, bty="l")
#terminando o save
dev.off()

					#105.2
#lendo arquivo
#esaligna<-read.table("esaligna.txt", head=TRUE, as.is=TRUE, sep="," )
esaligna

#criando a matrix pros boxplots do grafico 2
ligna<-matrix(nrow=6, ncol=2)
colnames(ligna)<-c("dap_mean","dap_sd")
rownames(ligna)<-c("16", "17", "18", "22", "23", "32")

#eh possivel indexar das duas formas seguintes:
ligna[,"dap_mean"]<-tapply(X=esaligna$dap, INDEX=esaligna$talhao, FUN=mean)
ligna[,2]<-tapply(X=esaligna$dap, INDEX=esaligna$talhao, FUN=sd)

#abrindo o arquivo pra salvar:
jpeg(filename = "Dispersao_105.2.jpg", width = 2000, height = 800, units = "px", pointsize = 12, quality = 100, bg = "white",  res = NA)
#primeiro plot
par(mfrow=c(1,2))
boxplot(dap~talhao, data=esaligna, xlab="Talhao", bty="L",ylab="DAP", main="Dispersao de DAP por Talhao")
text(x=3.501901, y=22.40602, "a")
#####

#segundo plot
par(pch=19)
plot(ligna[,1], xlim=c(1, 6), ylim=c(0,21),xaxt="n" , col="black", xlab="Talhao", ylab="DAP", main="Dispersao de DAP por talhao")
axis(side=1, at=1:6, labels=c(16, 17, 18, 22, 23, 32))
segments(x0=1, y0=(ligna[1,1]-ligna[1,2]), x1=1, y1=(ligna[1,1]+ligna[1,2]) )
points(x=1, y=(ligna[1,1]-ligna[1,2]), pch="-" )
points(x=1, y=(ligna[1,1]+ligna[1,2]), pch="-" )
segments(x0=2, y0=(ligna[2,1]-ligna[2,2]), x1=2, y1=(ligna[2,1]+ligna[2,2]) )
points(x=2, y=(ligna[2,1]-ligna[2,2]), pch="-" )
points(x=2, y=(ligna[2,1]+ligna[2,2]), pch="-" )
segments(x0=3, y0=(ligna[3,1]-ligna[3,2]), x1=3, y1=(ligna[3,1]+ligna[3,2]) )
points(x=3, y=(ligna[3,1]-ligna[3,2]), pch="-" )
points(x=3, y=(ligna[3,1]+ligna[3,2]), pch="-" )
segments(x0=4, y0=(ligna[4,1]-ligna[4,2]), x1=4, y1=(ligna[4,1]+ligna[4,2]) )
points(x=4, y=(ligna[4,1]-ligna[4,2]), pch="-" )
points(x=4, y=(ligna[4,1]+ligna[4,2]), pch="-" )
segments(x0=5, y0=(ligna[5,1]-ligna[5,2]), x1=5, y1=(ligna[5,1]+ligna[5,2]) )
points(x=5, y=(ligna[5,1]-ligna[5,2]), pch="-" )
points(x=5, y=(ligna[5,1]+ligna[5,2]), pch="-" )
segments(x0=6, y0=(ligna[6,1]-ligna[6,2]), x1=6, y1=(ligna[6,1]+ligna[6,2]) )
points(x=6, y=(ligna[6,1]-ligna[6,2]), pch="-" )
points(x=6, y=(ligna[6,1]+ligna[6,2]), pch="-" )
text(x=2.852165, y=20.05075, "b")
#####

#finalizando o save
dev.off()

					#105.3
#lendo os dados
arq<-read.table("arq_105.txt", header=TRUE, as.is=TRUE, sep=",")

#abrindo  o arquivo pra salvar:
jpeg(filename = "Dispersao_105.3.jpg", width = 2000, height = 800, units = "px", pointsize = 12, quality = 100, bg = "white",  res = NA)

#preparando o plot:
par(mfrow=c(1,2), lwd=2)

#primeiro plot:
par(bty="l", pch=17, family="serif")
plot(x2~x1, data=arq, xlim=c(0.5, 2.2),ylim=c(0,3), xlab="Log(Patch size)(ha)" ,ylab="Euclidean distances")
text(x=2.2, y=3,"a", cex=1.3)
segments(x0=0.499641, y0=2.189261, x1=2.147582, y1=1.168715)
#####

#segundo plot:
par(bty="l", tck=0.03, family="serif")
boxplot(x2~y2, data=arq, ylim=c(0,3.1), xlim=c(0, 7), outline=FALSE, xaxt="n")
points(x=c(1,2,4), y=c(3,3,3), pch=8 )
points(x=2.90, y=3, pch=8)
points(x=3.10, y=3, pch=8)
points(x=5, y=3, pch=8)
points(x=4.80, y=3, pch=8)
points(x=5.20, y=3, pch=8)
text(x=6.7, y=3,"b", cex=1.3)
axis(side=1, at=1:6, labels=c("Small", "Medium \nEdge", "Medium \nInterior", "Large \nEdge", "Large \nInterior", "Control"))
#####

#termianndo o save do plot
dev.off()

